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A two-dimensional small bias model has been developed for a patterned metal current 
collector | mixed oxygen ion and electronic conductor (MIEC) | patterned metal 
current collector electrochemical cell in a symmetric gas environment. Specifically, we 
compute the electrochemical potential distributions of oxygen vacancies and electrons 
in the bulk and near the surface for Pt|Sm .i5Ceo.850i.925|Pt symmetric cell in a H 2 — 
H2O — Ar (reducing) atmosphere from 500 to 650°C. Using a two-dimensional finite- 
element model, we show that two types of electronic current exist within the cell: an 
in-plane drift-diffusion current that flows between the gas | ceria chemical reaction 
site and the metal current collector, and a cross-plane current that flows between 
the two metal electrodes on the opposite side of the cell. By fitting the surface 
reaction constant k® to experimental electrode resistance values while fixing material 
properties such as bulk ionic and electronic equilibrium defect concentrations and 
mobilities, we are able to separate the electrode polarization into the surface reaction 
component and the in-plane electron drift-diffusion component. We show that for 
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mixed conductors with a low electronic conductivity (a function of oxygen partial 
pressure) or a high surface reaction rate constant, the in-plane electron drift-diffusion 
resistance can become rate-limiting in the electrode reaction. 

1 Introduction 

Mixed ionic and electronic conductors have received significant attentions for their 
potential as fuel cell components, permeation membranes, oxygen storage capacitors, 
electrochemical sensors, etc. Electrical d.c. and a.c. conductivity measurements 
of the bulk focused mostly on separating the ionic and electronic contributions to 
the electrical conductivity as well as determining the dielectric and chemical capaci- 
tances. Physically derived one- dimensional models have aided in the interpretation of 
electrical conductivity data in the bulk and materials defect chemistry. On the other 
hand, investigations of interfaces in mixed conductors, specifically the gas | mixed 
conductor and the metal | mixed conductor interface, require a two-dimensional 
model to adequately describe the system due to nonlinearity of the electrochemical 
potential lines near the interface. However, the majority of the work in the com- 
munity employs the one-dimensional treatment pQ, [2] and only a handful of works 
attempted to scale up to two dimensional models, [3] [I] and [5] . In this work, we de- 
veloped a two-dimensional small bias model for the patterned metal current collector 
| mixed conductor | patterned metal current collector cell. In particular, we focus on 
the numerical analysis of the cross-plane electronic current that flows through the 
mixed conductor between patterned metal stripes on both side of the cell, and the 
in-plane electronic current that flows between the gas | mixed conductor interface 
and the metal. In order for an electrochemical reaction to occur on the gas | mixed 
conductor interface, electrons need to diffuse from the reaction site to the external 
circuit and vice-versa. Such a step, termed in-plane electron drift-diffusion, could 
play a significant role in the interfacial behavior of mixed conductors, particularly 
for those exhibiting a low to moderate bulk electronic conductivity. 
For this study, we have selected Pt|Sm . 15 Ceo.850i.925(S , DC)|Pt as mixed oxygen 
ion and electron conductor model system. High oxygen ion conductivity of acceptor- 
doped ceria at intermediate temperatures (500 — 700°C) has attracted a great deal of 
interest in the solid oxide fuel cell (SOFC) community. In addition, under mildly re- 
ducing condition, doped ceria exhibits a moderate electronic conductivity (~ O.lS/cm 
at 650°C, po 2 — 10~ 25 atm [6]), making it attractive for fuel cell anode applications. 
Recent studies have also shown that, when operating SOFCs on hydrocarbon gases, 
ceria-based anode is significantly less susceptible to carbon coking [TJ. Insight into 
the in-plane electron diffusion path in ceria could lead to improved designs of anode 
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geometries and reduced interfacial resistance. 

2 Background 

The physical model, depicted in Fig. la, consists of a mixed oxygen ion and electron 
conductor (ceria) with patterned metal current collectors (Pt) on both sides placed 
in a uniform gas environment (H2 — H2O — Ar mixture). The patterned metal current 
collectors permit the system to be reduced to a repeating cell (Fig. lb) using mirror 
symmetry lines (r 1; r 2 , T 3 ). The thickness of the cell is given by 2Z 2 = 1mm, while 
the surface dimensions are 2W\ ~ 3^tm, the width of the metal | ceria interface (r 4 ), 
and 2W2 ~ 5/im, the width of the gas | ceria interface (T^) (Fig. 2). Two charge 
carriers species are considered: oxygen vacancies, denoted by the subscript ion, and 
electrons, denoted by eon. We solve the electrochemical potential and current of both 
charge carriers using a linear and time-independent model. 

We assume that the gas | ceria interface is the prevailing surface active site facilitating 
the reaction between electrons and oxygen vacancies in the oxide and the gas phase 
species. In other words, the gas | metal | ceria triple-phase boundary interface has 
a negligible contribution to the surface reaction. As shown later in this work, mixed 
conductivity allows electrochemical reactions to take place away from the triple-phase 
boundary. As a result, the metal | gas | ceria interface, a ID line, has substantially 
less area for reaction compared to the gas | ceria interface, a 2D area. We further 
treat the surface chemistry as one global reaction, and do not consider diffusion of 
adsorbed species on the surface. Combined with the final assumption that the metal 
I ceria interface is reversible to electrons, we are only considering two steps in the 
electrode reaction pathway: the surface reaction, and the electron drift-diffusion from 
the active site to the metal current collector. 

3 Model 

3.1 Governing Equations 

Mixed-valent oxides such as SDC exhibit mixed oxygen ionic and n-type electronic 
conductivity behavior under reducing conditions. Oxygen ion motion is attributed 
to vacancy hopping while electron motion is attributed to thermally activated small 
polaron hopping ([S] and [H]). This complex behavior of ceria at the atomic scale 
can be averaged out at the mesoscale (5nm and up) leading to constitutive equa- 
tions drawn from non-equilibrium thermodynamics [10] and [II]). For electrons and 
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oxygen vacancies we can write the following continuity equation: 



dt 



+ V-j m = fl m (1) 



where c m is the concentration of charge carriers m expressed in terms of number of 
particles per unit volume, j m is the particles mass flux, and Cl m is the net rate of 
generation, taken to be zero in our system since there is no source or sink of masses. 
The sum of all the charges in the system, — p = e (B — ^2 z m c m ) (where B is the 
concentration of the acceptor dopant), can be related to the electric potential, 0, via 
the Poisson's equation: 

V-(V</>) = -^ (2) 

We assumed that the permittivity, e, is a constant with respect to position. Finally, 
we can relate the mass flux to the carrier concentration and electric potential through 
the diffusion-drift equation: 

jm D m c rn ^ - — (3) 
k b T 

where fc& is the Boltzmann constant and T is the absolute temperature. In the dilute 
limit, the electrochemical potential fi m is given by: 

Am = aC + k bT log ( J + z m e<j) (4) 

where pP m is the standard chemical potential. In the steady state case and assuming 
that the diffusion coefficients are independent of position, combining Eq. [T] and [3] 
yield 



V ■ (c eon V/i eon ) = (5a) 

V • (CionVfiion) = (5b) 

and substituting jl m with Eq. H] while making the electroneutrality approximation 
(B + c eon - 2c ion = 0) further gives 

Ac eon - Vc eo „ • V0 - c eon A(j) = (6a) 

Ac eon + Vc eon ■ V0 + (B + c eon ) A0 = (6b) 

where the symbol A (•) indicates the operator V • (V(-)). 
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3.2 The Behavior of the Bulk 



We indicate the equilibrium quantities, such as electron and oxygen vacancy concen- 
tration, with the superscript (0). In order to determine equilibrium concentrations 
of charge carriers, we consider the following gas phase and bulk defect reactions: 

# 2 (gas) + 2 (gas) ^ # 2 0(gas) () 
O x Q ^ V£ + |0 2 (gas) + 2e' {) 

where the Kroger- Vink notation is used. We can also write the following equilibrium 
constants: 



K 9 = (8a) 
Ph 2 Po 2 

( (o) \ 2 (0) 

Kr = ( C -f) C ~fp£ (8b) 

(0) (0) 

1 + ^-2^ = (8c) 
B B 

Pk 

where pk = and pk is the partial pressure of species k. We solve the equilibrium 

latm 

concentrations of vacancies cfj n and electrons cion in the dilute limit at a given 
temperature and partial pressures. 

We suppose a small bias off equilibrium is performed in our system. Experimentally, 
this is achieved by subjecting the electrochemical cell to a small DC voltage relative 
to the open circuit voltage. Alternatively, we can obtain the same information by 
making an AC impedance measure at open circuit and taking the resistance at the 
the low frequency limit where the frequency approaches zero. We indicate these small 
perturbations with the subscript (1). These working conditions can be summarized 
as follows: 

• r - r (0) 4- r (1) with Vr (0) - and lr (0) I ^ lr (1) I 

• n - f/ (0) + f/ (1) with V/~/ (0) - and l// (0) I ^ I 

• peon — peon ~T~ fJ-eon Willi V \leon — U ailQ \jJ, e on\ \^eon\ 

• c ion = cfj n + with Vc££ = and \c^ n \ > \c\H\ 

• Vion = A + P-i! with V/i^ = alld l^ionl > 1/4™ I 
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+ 4> {1) with V0 (o) = and |0 (o) | > 



Substituting the above definitions and properties above into (jSJ) yield the following 
first-order relations: 

Au (1) - 

Au (1) = ( } 

We remark that linearization of the electrochemical potentials around equilibrium 
leads to: 



fl eon = k b T l0g HP - e ^ + /^eon 
V '-'eon / 

/2 ion = k b T log(^p) + 2e<p + fi° ion 

\ C ion J 



.(1) 



H'iori ~ " " ' (o) 



%r (io) 



ion 



Furthermore, applying the electroneutrality approximation to the equilibrium as well 
as the first-order perterbation in carrier concentrations gives 



R _l_ „(°) 

£1 = (Ha) 



(1) Ceon 



(1) 

(lib) 



The [IT] allows us to reduce Eq. [10] to a non-singular linear transformation between 
the variables I <p^\ \ and I fish, /4on r- This implies we could use either set of 
unknowns to fully describe the problem. In the reminder of the paper we will use 
the following two quantities rS^ = and (fi^ = ^ttf-- 

C-eon " 

3.3 The Behavior of the Boundary 

Realistic boundary conditions are complex due to the formation of charge double 
layers [12], [13]. Work by Fleig et al. suggests that such electrification effects is 
relevant for SOFC mixed conducting cathodes [H], [15]. For simplicity, we do not 
consider charge double layer in our system. Furthermore, we do not consider surface 
diffusion as the the need to specify the surface roughness may lead to over-fitting of 
the data. 
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As shown in Fig. 1, there are five boundaries in our electrochemical system. We start 
with the simplest boundary conditions. It follows from symmetry that d x jx eo \ = 0, 
®xfiion = on T 2 and T 3 . Since the metal is ion blocking we will have dyjlf^ = on 
T 4 . By assuming that response of the metal to any perturbation is fast compared 
to the oxide, we can take (ffi*' uniform on T 4 . Thank to linearity we can choose 
0« = k b T/e on T 4 (so that = 1) and = on T x . 

The remaining boundary, r 5 , is complex due to the gas-solid surface reaction electro- 
chemistry. Specifically, the fuel cell anode condition under which our computation is 
performed requires us to consider the interaction of oxygen vacancies and electrons 
SDC and gas-phase hydrogen, oxygen and water vapor. A complete treatment of the 
surface require a detailed understanding of the electrochemical reaction pathway and 
kinetics parameters of various reactions that occur in series and/or parallel. How- 
ever, there is little experimental data in literature regarding the surface reactions 
for SDC (or for any other composition of doped ceria). In the case of ceria, AC 
impedance spectroscopy is unable to separate multiple processes that occur on the 
electrode | electrolyte interface due to the overwhelming " chemical" capacitance that 
results from redox of cerium cation between +4 and +3 oxidation states. 
In this work, we treat the surface reaction pathway as a single reaction. Specifically, 
we assume that the surface chemistry can be described by: 

# 2 (gas) ^ H 2 (gas) + V*' + 2e' (12) 

Furthermore, we assume that the rate of reaction, specifically, the rate of injection 
of vacancies at T 5 satisfy 

lion ■ e y = k f p H2 - k r p H2 oc ion c 2 eon (13) 



where e y is the unit vector that is perpendicular to Ts, kf and k r are the forward and 
reverse reaction rate constants, respectively. At equilibrium, the net rate of injection 
of both oxygen vacancies and electrons are zero, so kf and k r can be related to the 
equilibrium concentrations of the reactants and products: 

kr = (14 ) 



Ceon + D I I Ceon I PH 2 



Using the same approach as Section 3.2, we compute the perturbation in the bound- 
ary condition upon applying a small bias perturbation. Combining [13] and [TH and 
perturbing c eon = c^Jn + c£Jn, we obtain: 
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j^-e, = 2k f 1 + ^ U Ha „« (15) 



We go a little further and suppose kf — 2 '°" fcy fc/ = ^/Po 2 x 4L H — ' w ^ ere we 

choose 8 = —1/4. 



It is important to note that the choice of 8 is based on the oxygen partial pressure 
dependence of the rate limiting step(s) in the surface reaction. Since identifying the 
rate-limiting step in the surface reaction is beyond the scope of this work (as we 
described the surface reaction with a global reaction), we selected the 8 value so 
that the po 2 dependence matches the experimental results [6] that will be used for 
data analysis. 



3.4 Formalization of the Model 



If one defines x = Z c x (l c = lOfim) ,A^, = —k 



f ,(o) 



1 - 



D; 



and A r 



^ion 



D c (0) 

j _|_ ^_ion^ion 



,(°) 



eon^eon 



then ([9]) and the boundary conditions of the previous section can be summarized as 
follows: 

A^ 1 ) = 
A^nW = 



(16) 



f 0« = 

<9*0 (1) = 
0« = 1 



A = 
A dsn^ = 



(0) 



(0) 

'■ -co n 



{ dy<j>V> = A^ A djjnW = A n n« 



* The units of the kf's 



[kf] = 



^particles 

s x m 2 
^particles 



on r\ 

on r 2 

on r 4 

on r 5 



a r. 



(17) 



let us look at order of magnitude of kf. po 2 = 10 24 , h = 10 5 m, D ion = 10 10 m 2 /s and fe9 w 10 32 , 
so kf w 10 



32 w 10 



=r x 10" 



=10 



21 # particles 



10 -3m^_ 10" 7 ^4 

— cm 2 
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3.4.1 Numerical Method 



In order to solve numerically the Equations (|T6l) with the boundary conditions ( fTTl) 
we employ an h-adapted finite element method (FEM). FEM is well known for elliptic 
problem such as the one we are studying, see for example [16] and [T7] . 
In order to employ FEM appropriately, we first recast the problem in the following 
weak form, where m and ip are test functions on the domain Q: 

I Vm • Vn (1) d^ + 4^gL ! m<%0 (1) d 7 - A n j mn w d 7 

JCl Ceon J T4 JVs 

j Vi) ■ V0 (1) dA-A^ I ijn {1) d 7 

with the addition of the conditions that: 

f <p) = A = 1 on r x 

\ = 1 on r 4 

The discrete version of equations (|18p are then solved using FreeFem++ [18]. The 
equations are initially discretized on a triangular unstructured mesh, using quadratic 
continuous basis functions with a third order bubble. The mesh is adaptively refined 
up to seven times at each solution step and the a posteriori adaptation is performed 
against jlem- The h-adaptation ensures high regularity of the H 1 a posteriori esti- 
mator [19], locally below 0.01%, and it guarantees that the mesh is finer where the 
sharpest gradients occur. We note that mesh adaptivity results is coarseness every- 
where except in the vicinity of the interfaces, in particular the refinement increases 
as we approach the triple phase boundary; this fact indicates strong nonlinearities 
around that area. Eleven integral tests were also implemented in order to ensure that 
at solution step the numerical method is consistent with the boundary conditions and 
it globally satisfies conservation of charge. Finally we note that FreeFem+-|- execu- 
tion time is comparable to custom-written C++ code and its speed his enhanced by 
the utilization of fast direct linear solvers such as the multi-frontal package UMF- 
PACK [20] • Due to the sparsity of the problem we make extensive use of this last 
feature. 

3.5 Value of the Polarization Resistance 

To compute the electrode polarization resistance, let us first consider the relevant 
electrochemical currents that exist within our system. Due to the mixed conducting 
nature of ceria, there will be an inherent cross-plane electronic current, termed l£ p , 



= 

(18) 

= 
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that flows between the metal current collectors located on the opposite side of ceria 
(Figure [2]). The surface reaction taking place on Ts will simultaneously inject one 
oxygen vacancy and two electrons into ceria. Two distinct current result: the cross- 
plane ionic current Ig F that flows between the two sides of the cell, and the in-plane 
electron drift-diffusion current Ig P that flows between gas | ceria interface and the 
metal | ceria interface. By electroneutrality, = Ig P . We formally define these 
currents as follows: 



Ie P = / Jeon • n dx = / j eon • n dx 



(20) 



Ig P = / Jeon • n dx = / j eon • n dx 

Note that to obtain I IP , we could integrate the current density either over T 5 or 
over <^f(r 5 ), which, as depicted Fig. El represents some fraction of T 4 accessed by 
the current injected from the gas | ceria interface. Direct comparison with the work 
of Jamnik and Maier, p], and Lai and Haile, [6], leads to the following definitions 
of the bulk electronic resistance R eon , bulk ionic resistance Ri on , and the electrode 
polarization resistance normalized by the sample area R^ on : 



Reon 2 



< f^eon >r 4 — < ^eon >Ti < Peon >r 4 



^ _ Q < ^eon > r 5 ~ < Veon ^£1 _ ^ < Heon >r 5 
^ 9 „ , J 'g 



R; 



± . < ^eon >r 4 _ < f^ion >r 5 



lull 



3g 



jCP JIP 

where j e = Wl e +W2 and j g = Wl 9 +W2 are the current densities averaged over the total 
sample area. Our two-dimensional model allows us to separate various contributions 
of Rj~ on - We shall discuss two physically significant separations of the contributions. 
First, it is possible to separate R^ n into a surface reaction resistance term R sur f and 
in-plane electron drift-diffusion term R eon -DD- Rsurf corresponds to the electron- 
vacancy electrochemical potential difference at the gas | ceria interface, and R e0 n-DD 
corresponds to the electron potential difference at the gas | ceria interface and at the 
metal | ceria interface. 

d _ < Kon >r 5 ~ < Km, >r 5 _ A , n \ < n^> > Fs 
SUrf ~ ^ h "V W j 9 (22) 

Reon—DD Rion Reon— ion 
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Physically speaking, R sur f represents the resistance associated with the chemical 
transformation of electrons to vacancy at the gas | ceria chemical reaction site. Specif- 
ically, Rsurf is due to the migration of oxygen vacancy from the bulk to the surface, 
and the subsequent chemical reactions that give rise to the electronic current. Fi- 
nally, the drift-diffusion of the injected electrons from the reaction site to the metal 
current collector results in R e0 n-DD- 

Alternatively, we could also separate R^ n into a "true" polarization term, R po i, and 
a "deviation" term, R avg , that results from averaging: 

p ^ fieon ^V 3 ! ( r 5) _ fiion 

-LLpol 



h (23) 
Ravg *~^ion *~^pol 

In Eq. (121]) . R^ n is proportional to the difference of the electronic electrochemical 
potential averaged over the metal | ceria and averaged over gas | ceria interface. In Eq. 
( |23|) . we define the R po i by averaging only some portion of (rather than over the 
entire interface) by considering the interface mapped by current lines injected from 
gas | ceria interface. R aV g, defined as the difference between the electrode polarization 
resistance and the true electrode polarization, is simply a spurious contribution due 
to averaging. 



4 Results 

4.1 Potential Distributions and Surface Regions 

Electrochemical equipotential lines for oxygen vacancies (Fig. El right) calculated 
using various values for the surface reaction rate constant, k®, reveal that the poten- 
tial and current distribution exhibit a relatively weak dependence on k®. In general, 
oxygen vacancy equipotential lines bend as they approach the oxygen vacancy block- 
ing metal|ceria interface (r 4 ) from the bulk. On the other hand, equipotential lines 
for electrons (Fig. El left), display substantial deviations from those for oxygen va- 
cancies, due to the presence of two current sources: cross-plane electronic current 
that flows between the current collectors on opposite side of the cell, and the in-plane 
electronic current injected by the surface reaction that flows between the metal|ceria 
(r 4 ) and the gas|ceria (r 5 ) interface. The electron potential distributions also depend 
strongly on the magnitude of kp indicating that electronic current injected from the 
surface reaction taking place at (Fa) strongly influence the electron penetration depth 
of the so-called "surface region." 
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Fig. [6] shows the boundary of the surface region, given by the "trajectory" of electrons 
injected from the surface reaction site furthest from the metal current collector (the 
intersection of and r§ in Fig. [1]). Physically, the surface zone can be viewed as 
a region where electronic current is entirely the in-plane electronic current {Ig P ), 
rather than the cross-plane current I^ p '. The surface region dimensions (Fig.s [31 El 
[7j) are specified by the largest length (/), largest depth (d), and the area (A). All 
dimensions increases as a function of k®. As the penetration area increases, the in- 
plane electrons will flow through a larger cross-section of ceria, thereby reducing the 
diffusion resistance. It is interesting to note that the surface region approaches an 
asymptote for large k®, suggesting that when surface reactions are sufficiently fast, 
i.e. when they are in electrochemical equilibrium, the total electron injection current 
will be dominated by the in-plane electron diffusion resistance. 
The surface region dimensions also grow with increasing po 2 , though it is more pro- 
nounced for higher k®. At lower values, the penetration area is virtually indepen- 
dent of po 2 ■ The penetration depth is a function of the relative magnitude of I\ p to 
. As the ratio I IP /I^ p grows, for instance, as a function of po 2 , the penetration 
depth is expected to increase. In Fig. d we see that an increase in the penetration 
dimensions is indeed accompanied by an increase in I IP /I^ P . 



4.2 Electrode Polarization Resistance 



Eq. ( 1231) states that the electrode polarization resistance, Rj~ on , can be expressed 
as a sum of the true polarization term R po i and a deviation term R aV g that results 
from averaging the electrochemical potential of electrons across the entire metal | 
ceria interface T 4 rather than just the region accessed by the in-plane electronic 
current (pf{T 5 ) (Fig. [3]). We examine the extent of deviation of the R^ n from R poi 
by computing: 

fpol = I? (24) 

ion 

(25) 

Under a variety of conditions, f po \ is very close to unity (Fig. [10] top), indicating that 
the deviation term is quite small compared to the true polarization resistance. For 
the remainder of our analysis, we approximate R po i = R^ n - 



4.3 Electron Diffusion Resistance 

In Eq. (|22|) . the total electrode polarization resistance, R^ n , is expressed as a sum of a 
surface reaction resistance term, R sur f, and an in-plane electron diffusion resistance 
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term, R eon ^£>D- Under the moderately reducing po 2 regime (where the electron 
carrier concentration is negligible compared to the extrinsic oxygen vacancies formed 
by acceptor-doping), R eon _DD is proportional to approximately Pq^ (Fig. [9] bottom), 
following the same po 2 dependence as bulk electronic resistivity As for R sur f, it is 
also proportional to p 1 ^ as a result of our choice of (3 (Fig. [9] top). Turning to the 
dependence on k®, we observe that both R sur f and R eon -DD decreases with increasing 
kj. However, a significant difference between R eon _DD and R sur f is that, in a log-log 
plot the former approaches an asymptotic value as a function of k®, whereas the 
latter does not. This interesting behavior of R eon -DD is directly related with the 
asymptotic behavior of the the penetration depth of electrons injected from r§ to T4 
(and vice- versa) as a function of k® (Fig. [6]). To help us examine Ig P , we further 
define fractional surface reaction resistance and fractional electron drift-diffusion as: 

£ ^surf 

Jsurf ~ (26) 

feon—DD 1 fsurf 

Plotting f sur f as a function of k'j (Fig. [10] bottom) reveals that when the surface 
reaction is very fast (i.e. large k'j) , f sur f approaches zero and R^ n is dominated by 
Reon-DD- On the other hand, when the surface reaction is slow, f sur f approaches 
unity and i?^ n is dominated by R sur f, as would be expected. When considering only 
material property dependencies (i.e. neglecting po 2 , T, and sample geometry) , R sur f 
is only a function of kj whereas R eon -DD is a function of both kj and a eon . As kj 
tends toward infinity, R sur f approaches zero and R eon -DD approaches an asymptotic 
limit that is a function of only a eon . In other words, as the surface reaction resistance 
term becomes negligible, electron carrier concentration and mobility alone determines 
the penetration dimensions and therefore R^ on - The condition under which f sur f 
approaches zero corresponds to the physical case where the electrode reaction is 
limited by the rate in which the electrons migrate from the gas | ceria reaction site 
to the metal rather than the rate of surface reaction. Generally speaking, for a 
wide-bandgap mixed conductor exhibiting a low or moderate electronic conductivity 
and high k'j, such as ceria, in-plane electron drift-diffusion cannot be neglected. 
Accordingly, the electron diffusion length, (separation between the metal in Fig. [T] 
top), needs to be tuned in order to minimize the electrode polarization resistance. 

4.4 Topological Considerations 

There are two degrees of freedom in the metal current collector topology: the metal 
stripe width (2Wi) and the intermetal distance (2W 2 ). Fig. [TT] show parametric plots 
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of the the fractional surface resistance (top row) and the total electrode polarization 
resistance (normalized for the total sample area) (middle row) and as a function 

W 2 

of Wi,— — , and k f at select temperatures and po 2 - We observe the general trend 
W 2 

that increasing — — (gas | ceria interface to metal | ceria interface ratio) leads to 
W\ 

a reduction in the polarization resistance. Specifically, under the conditions that 

the fractional surface reaction resistance is greater than 0.99, we observe a linear 

W 2 

decrease in the polarization resistance with increasing — — . When the electron drift- 

W\ 

diffusion resistance is negligible, the current density of electrons injected from T 5 

is essentially uniform as a function of the position (and distance to the metal), 

and thus the fraction of area available for electrochemical surface reaction, given 
W 2 

by / = ^ — , determines the polarization resistance. Fig. [TT] (bottom row) 

shows the polarization resistance normalized by the gas | ceria interface area and 
confirms that the normalized resistance remains relatively constant as long as the 
electron diffusion fractional resistance is negligible. However, as the fractional surface 
resistance decreases (due to an increase in k% for instance), electrochemical surface 
reactions taking place closer to the metal will inject a larger current into the oxide. 
For example, at = 7.5 x 10 33 , po 2 = 4.1 x 10~ 26 and T = 650°C, a significant 
nonlinearity as well as a distinct minima in the polarization resistance as a function 
of W\ and W 2 (Fig. [TTJis observed (top right)). Furthermore, Fig. [TT] (lower right) 
shows the polarization resistance normalized for the gas | ceria interface area begins 
to deviate from the constant values, confirming parts of the interface is becoming 
less active due to increased electron diffusion resistance at spatial positions further 
away from the metal current collector. 

In general, when the surface reaction rate constant is small or when the bulk elec- 
tronic conductivity is large, one should increase / in order to increase maximize the 
area available for surface reactions, as long as the electron diffusion fractional resis- 
tance is kept low. On the other hand, when the surface rate rate constant is large or 
when the bulk electronic conductivity is small, one needs to find an intermediate / 
in order to balance the area available for surface reaction and the in-plane electron 
diffusion distance. 

4.5 Comparison to Experimental Results 

We fit the polarization resistance data obtained by Lai and Haile [6] using AC 
impedance spectroscopy on a cell geometry consistent with our model description. 
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The experimental result was based on a porous Pt\SDC\ porous Pt cell in H2 — H2O — 
Ar. We approximated the porous Pt electrode as line patterns by estimating W\ and 
W2 based on the actual pore size and interpore distance. We fit the polarization 
resistance using as the only parameter and fixed dopant and equilibrium car- 
rier concentration according to the values obtained in the experiment. It should 
be noted that all parameters were obtained from the same electrochemical cell by 
Lai and Haile and are highly self-consistent. The fitting (Fig. IT2T) shows computed 
corresponding to the polarization resistance obtained experimentally. Because we 
phenomenologically set thepo 2 dependence of R sur f to 1/4 so that Rj^ n would exhibit 
the same po 2 as the experimental data (and the other component of the polarization 
resistance, R eon -DD, is also proportional to pQ^ in the same way as the bulk elec- 
tronic conductivity), obtaining the same dependence in po 2 for the experimental and 
fitted value is automatic. 

Taking the fitted k'j values, we can further separate the polarization resistance into 
the surface reaction and the electron drift- diffusion contributions. At the tem- 
peratures and po 2 examined, the computed f sur f (Fig. [Hj) is close to unity (for 
W\ ~ 1.5/xm and Wi/W\ ~ 1-67), implying that the surface reaction step is the 
rate-limiting step. To examine the dependence of f sur f on the geometric parameter 
(which directly influences the electron diffusion length and area of the gas | ceria 
interface), we fit k® to the polarization resistances while varying W\ and W2 ■ The 
parametric plot (Fig. [TJJ again shows that f SU rf is close to unity for a wide range of 
W±, W 2 , T and po 2 - However, we do observe the general trend that f sur f decreases 
slightly with increasing W\ and decreasing W 2 /Wi. Decreasing W%jW\ (at a fixed 
W\) reduces the electron diffusion length and reduces the area of the gas | ceria 
interface, and fcy needs to be increased in order to fit to the observed polarization 
resistance (Fig. [13]). For the same reason discussed in Section 4.3, this leads to a 
decrease in f sur f. On the other hand, increasing W\ (at a fixed W 2 jW\) increases the 
electron diffusion length without affecting the available reaction area. As a result, 
increased R eon -DD leads to an decrease in f sur f. 

Approximating a grid-like porous metal on ceria as line patterns could lead to some 
errors, such as over-estimating the fraction of gas | ceria interface and the electron 
diffusion length. However, given that the computed f sur f is far from 0.5 (the case 
where surface reaction and electron drift-diffusion are equally co-limiting) for a wide 
range of W\ and W 2 , these errors will not change the f sur f significantly and will only 
re-scale the magnitude of the resistances slightly. Therefore, based on the numerical 
analysis in this work, the electrode reaction in porous Pt^DCI porous Pt cell in 
H 2 — H 2 — Ar is likely to be surface reaction limited. 

Finally, it should be noted that our assumption that the electron mobility and equi- 
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librium carrier concentration is the same in the near-surface region and in the bulk 
directly determines the contribution of the in-plane electron drift-diffusion resis- 
tance to the electrode polarization resistance. Since electron penetration depth is 
predicted to be on the order of 1/xm, dopant segregation and presence of blocking 
grain boundaries near the surface could in principle affect the local electron mobility 
and concentration. 

5 Conclusions 

A two-dimensional electrochemical model has been developed for mixed conductors 
with patterned metal current collectors. Numerical simulation for a Pt | SDC | Pt 
in reducing atmosphere revealed a strong nonlinearity in the electronic potential 
and current distributions near the surface. In particular, we show that the in-plane 
electron drift-diffusion current plays a crucial role in determining the surface elec- 
trochemical behavior. Under certain conditions, the in-plane electron drift-diffusion 
resistance could dominate the electrode resistance. 
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A Definitions 

Some nomenclature: 

1- c e on = concentration of electrons in #particles/m 3 

2- c ion = concentration of vacancies in #particles/m 3 

3. B = concentration of negatively charged background particles in #particles/m 3 

4. /i e0 n = chemical potential of electrons = log ( %p- ) + /ig 0n 

\ c eon J 

5- fJ>ion = chemical potential of vacancies = log ( %^ ) + n° ion 



6- Jieon = electrochemical potential of electrons = /i eon — e0 

7. fiion = electrochemical potential of vacancies = /ii on + 2ecf> 

8. jl* on = star electrochemical potential of electrons = — -/i eo 

9. Ji\ on = star electrochemical potential of vacancies = + ^//j, 
10. = electric potential 

11- ieon =electron current = eD eon c eon V ^ 

12. iion = vacancy current = -2eD eon c eon Vj^ 




k B T 



13. j 




eon 




vacancy ffux= —D ion Ci, 



■ion 



k B T 



15. U T = 



e 
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B Assumptions 

For our Pt | SDC | Pt cell, we make the following assumptions: 

1. the conditions are steady-state, i.e. d t (-) = for all unknowns; 

2. electroneutrality is satisfied (i.e. c eon — 2ci on +B(x.) 0) throughout the sample; 

3. the drift-diffusion equations describe correctly the fluxes; 

4. diffusivities of electrons and vacancies are constant and uniform; 

5. surface diffusion is frozen; 

6. the chemistry is correctly described in one global step. 
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Table 1: Insert table caption here 



T 


500°C 


550°C 


600°C 


650°C 




m? 
V 2 s 


4.762 x KT 8 


6.257 x 10" 8 


6.873 x 10~ 8 


8.123 x 10~ 8 






m 
V 2 s 


1.166 x 1(T 9 


2.070 x 10~ 9 


3.359 x 10- 9 


4.936 x 10~ 9 



K 9 

K r 



5.059 x 10 27 4.814 x 10 25 
5.008 x 10~ 22 2.263 x 10~ 20 



7.757 x 10 23 
6.610 x 10~ 19 



1.944 x 10 22 
1.340 x 10~ 17 



Wi \pm] 


w 2 H 


L [jum] 


£[21J[22_ 




" #particlcs 




1.5 


2.5 


500 


3.47 x 10 27 



21 



Fig. 1 Caption of Figure 1 
Fig. 2 Caption of Figure 2 
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Figure 1: top: A symmetric cell with patterned Pt stripes on both sides of dense ceria 
placed in a uniform gas environment, bottom: Schematic depiction of the boundaries. 
Ti, Y21 and 1^ are symmetry lines, while is the metal | ceria interface, and is 
the gas I ceria interface. 2W\ is the width of the metal, 2W-2 is the width of the ceria 
directly exposed to the gas phase, and 2Z 2 is the thickness of ceria. 
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Figure 2: Various electronic and ionic current within the cell. Solid line indicates 
the electronic current and dashed line indicates the oxygen vacancy (ionic) current. 
The superscript IP refers to the in-plane current that flows between the gas | ceria 
chemical reaction site and the metal current collector, while CP refers to the cross- 
plane current that flows between the metal current collectors located on the opposite 
side of ceria. 
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BULK REGION 



Figure 3: Illustration of the surface region, where the in-plane electronic drift- 
diffusion current prevails. The dimension of the region is indicated by the length 
I anddepth d. (pf(T 5 ) depicts the fraction of the metal | ceria interface mapped by 
the electronic current injected from gas | ceria interface. 
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Figure 4: Simplified one-dimensional equivalents circuit employed by Lai and Haile, 
and Jamnik and Maier, where R eon is the bulk electronic resistance, Ri on is the bulk 
ionic (oxygen vacancy) resistance, and R^ n is the electrode polarization resistance 
normalized by the cell area. /i* on and /i* on are the electrochemical potential of oxygen 
vacancies and electrons, respectively. 
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Figure 5: Electrochemical equipotential lines (left) and the corresponding current 
flow lines (right) computed for various surface reaction rate constants kfO at 650C, 
p02 = 4.1E-26 atm. 
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log 10 p O2 = -25.39, T = 650°C 
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Figure 6: The boundary of the surface region, where electrons undergo in-plane drift- 
diffusion between the gas | ceria inteface and the metal current collectors, computed 
for various surface reaction rate constants kfO at 650C, p02 = 4.1E-26 atm (top) 
and 2.1E-21 atm (bottom). 
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Figure 7: Dimensions of the surface region as a function of po 2 and at 650°C. 
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Figure 9: Absolute value of the surface reaction resistance R sur f (top) and the elec- 
tron drift-diffusion resistance R eon -DD as a function of p02 and ~kf at 650°C. 
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Figure 10: Fractional true polarization resistance (top) and fractional surface reac- 
tion resistance (bottom) as a function of k^p at 650°C, parametrized with respect to 
iogioPcv 
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Figure 11: Parametric plots of the fractional surface reaction resistance (top row), the 
electrode polarization resistance i?^ n normalized for the sum of the metal | ceria and 
gas | ceria interfacial area (middle row), and normalized for the gas | ceria interfacial 
area (bottom row), as a function of W2/W1 and W\ 
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Figure 12: Electrode polarization resistance , R^ n , plotted as a function of po 2 and 
A;^ at 500°C, 550°C, 600°C and 650°C. The open triangles show the experimental 
data obtained in ref 
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Figure 13: If we assume that that our initial choice of W\ and W 2 is not the correct 
one, it is legitimate to ask the following question: "which k® fit the ASRP data 
best?". We find that the fitting depends only on the ratio and not on the chosen 
value of Wi, the label indicates that W± = 0.5, 1.0, 1.5, . . . 3.0/im; subsequent pictures 
show that the Area affected by the surface T 5 goes with Wf. This is a check of how 
good the numerical study is. 
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Figure 14: Fractional surface reaction resistance, obtained after fitting kj to the 

experimental data in ref ???, and plotted function of p02 kf at 500°C, 550°C, 
600°C and 650°C 
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